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Abstract 

We study the effect of infrared renormalons upon shape variables that are 
commonly used to determine the strong coupling constant in e~^e~ annihilation 
into hadronic jets. We consider the model of QCD in the limit of large nj. 
We find a wide variety of different behaviours of shape variables with respect 
to power suppressed effects induced by infrared renormalons. In particular, 
we find that oblateness is affected by 1/Q non-perturbative effects even away 
from the two jet region, and the energy-energy correlation is affected by l/Q 
non-perturbative effects for all values of the angle. On the contrary, variables 
like thrust, the c parameter, the heavy jet mass, and others, do not develop 
any 1/Q correction away from the two jet region at the leading nj level. We 
argue that 1/Q corrections will eventually arise at subleading nj level, but that 
they could maintain an extra as{Q) suppression. We conjecture therefore that 
the leading power correction to shape variables will have in general the form 
C(g{Q)/Q, and it may therefore be possible to classify shape variables according 
to the value of n. 
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1. Introduction 

Tests of QCD carried out at e'^e~ colliders have received a considerable boost at 
LEP and LHC where, because of the large centre of mass energy, the pertur- 
bative character of jet production becomes quite prevalent. Although jet studies at 
LEP provide convincing evidence of the validity of the perturbative approach, the 
determination of the strong coupling from jets cannot be considered as solid as other 
type of determinations, like those from the Z hadronic width or deep inelastic scat- 
tering experiments [0]. In fact, even at LEP energies, there are substantial power 
suppressed effects that are corrected for using Monte Carlo models. The estimate of 
the theoretical error associated with these corrections is very difficult, and inherently 
model-dependent. It would be very desirable to acquire some knowledge of power 
corrections from theory alone, without the need to resort to models. 

Some sources of power suppressed effects are in fact understood as originating from 
factorial growth of the coefficients of the perturbative expansion arising either from 
the large momentum region (UV renormalons) or from the low momentum region 
(IR renormalons) of a certain class of feynman graphs (see ref. and references 
therein). IR renormalons, UV renormalons and instantons are the only known sources 
of factorial growth in the perturbative expansion. Instantons are known to give 
corrections that are suppressed by a very high power of the hard scale involved in 
the process, while renormalons may give 1/Q corrections to certain quantities. It 
has been argued in refs. |^ and |^ that 1/Q power corrections arising from infrared 
renormalons are present in certain jet shape variables, and that these corrections may 
also be described in a common framework in terms of a "frozen" running coupling 
constant. In ref. the strongest suggestion is made that the 1/Q corrections may 
factorize, and it may even be possible to describe power corrections in Drell-Yan pair 
production and in jet events in a unified framework. Power corrections to jet shape 
variables were also considered in a simplified model in ref. 0]. 

In the present work, we actually compute the effect of renormalons on jet shape 
variables in QCD in the limit of large rtf In this limit the theory is not asymptot- 
ically free, but we will try to infer the properties of the full theory just by changing 
the sign of the first coefficient of the beta function at the end of our calculation. Our 
attitude is that QCD is at least as bad as this limit. 

The remainder of this paper proceeds as follows. In Section 2 we will give an 
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introductory description of our calculation, without going into any technical detail. 
In fact, the essence of the physical picture that we develop is contained in this section. 
In section 3 we give a description of the full calculation, and deal with the subtleties 
associated with canceling the real and virtual divergences by performing the calcula- 
tion in dimensional regularization. In Section 4 we discuss the term (which we call 
"Sudakov" term) that has a factorized 3-body form, which we identify with the term 
discussed in refs. [H and In Section 5 we discuss the non-factorizable piece. We 
show explicitly that this term is different for thrust and for the heavy jet mass, which 
are quantities that have the same definitions at the 3-parton level, but differ at the 
4-parton level. In Section 6 we discuss the possible effects of subleading corrections, 
and what could be expected in the full QCD theory. Finally, in Section 7 we give our 
conclusions. 

2. Infrared renormalon in the large Uf limit 

We will first examine the effect of infrared renormalons in the limit of large nf. 
In this case the dominant graphs are those given in fig. |l|. Notice that there are two 



Figure 1: Dominant diagrams for e~^e into jets in the large Uf limit. 

types of contribution, one with three partons in the final state and one with four. 
When taking the square of the amplitude with four final state partons the fermions 
coming from the gluon splitting should not interfere with those coming from the 
photon vertex in order to give a dominant term in the large nj limit. 
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We will now discuss what we expect for the result of the calculation of jet shape 
variables. The aim of this discussion is simply to give the flavour of how the exact 
calculation, which will be presented in the following section, works. All details of the 
infrared cancelation will be dealt with in the next section. 

In some appropriate renormalization scheme the inclusion of all vacuum bubbles 
will amount to the replacement 

where 

and Tr = nf/2. The contribution of the graphs with three partons in the final state 
will be simply given by the Born cross section, with the replacement of eq. (|2.1| ) 
at A;^ = 0. In this limit our expression vanishes, which is to say, we have no virtual 
graphs after we have resummed the whole perturbative expansion. It does not vanish, 
however, order by order in perturbation theory, but has instead an expansion with 
infrared divergent coefficients. In the next section we will show that these divergences 
are canceled order by order in perturbation theory by the real term. For the sake 
of the present illustrative discussion, we will instead accept the fact that it vanishes, 
and concentrate on the finite remainder coming from the real process. 

For the corresponding real process, the amplitude will in general have the form 



1 + asbo (log ^ + 



ITT 



(2.3) 



where represent here the whole of phase space, except for the virtuality of the 
gluon, k^. In this expression we have factored out the infrared divergent term dk^/k^, 
so that A{k^, 0) is in fact regular for k^^O. Observe that since k'^ is positive, the 
logarithm in the running coupling acquires an imaginary part. Let us now suppose 
that we are computing some infrared safe shape variable S{k'^,(j)). Infrared safety 
implies that in the limit of small k'^ S{k'^,(j)) goes continuously to its three-body 
form, so that the cancelation between real and virtual infrared divergences takes 
place. We define 

G{k^)= f A{e,(j)) S{k\(j)) d(j). (2.4) 
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The value of 5* will be given by 



S 



A;2 



die 



I + as bo (log g + ztt) 



(2.5) 



which we rewrite as 



S 



az 



For the integral in the second term we get 



(2.6) 



G'(O) 



de 



at 



G(0) 



arctan 



"l + aA logg^ 



arctan 



1 + agfeo log 



GjO) 



arctan 



asboTT 



1 + as6o log % 



arctan 



asboTi 
asboTT 



1 + asbo log 4 



(2.7) 



where A is an infrared cutoff. There is a subtlety in the second step, where we use 
the identity arctan(x) = ±7r/2 — arctan(l/x), in which the 7r/2 takes the same sign 
as X. Thus our substitution is only correct when the arguments of both arctangents 
have the same sign. This condition is violated if bo is positive, and it might appear 
that we have neglected a term proportional to vr, with no coupling in front. However, 
we should remember that we are computing a perturbative expansion, and that our 
algebraic manipulations should always be interpreted as an order by order expansion 
in as- We should therefore always be reasoning by assuming that terms with factors 
of as are small, even if multiplied by infrared or ultraviolet divergent coefficients. In 
this sense eq. ( p.7|) is correct for either sign of bo- 

The infrared divergent term we obtain cancels against the virtual diagram. The 
cancelation is explicitly shown in the next section. Here we just assume that it will 
take place. In the present context, the real infrared divergent term vanishes in the 
same sense in which the virtual term was vanishing. We see therefore that the only 
place where we can obtain an infrared renormalon is the first term of eq. (^^.6]). Let 
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us assume that G{k'^) — G{0) oc for small k, and consider the integral 



dk'^ ( k 



a. 
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Jo 
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boTT 
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boll 



(l + asbo log ^ - asboz^ + (asboT^y 



/ dze-P'/^- 
Jo 1 



etc 



+ as6o log - asboz - iagbon 



ZQ+OO 



^^'g-p(.'-.0)/2 



as 



zo 



1 - asboz' 



(2i 



where 2; = log ji^, z' = z + zq and 2:0 = ivr — log This becomes 



fc2 



1 



1 



Im 
Im 



asboz' 

as 



agboz' — ie 



(2.9) 



and the first term is analytic in as- Rescaling z' we finally get 



boTT 
1 



Im 
Im 



ipn/2 ( 



ac 



1 - asboz' 
^ ^^exp{-z/as) 



p 

2bo 



z — te 



(2.10) 



The first term is analytic, while the second term has an infrared renormalon located at 
p/(26o), which corresponds to a power correction of the order of l/Q^. Observe that 
for positive bo we have found a definite prescription to bypass the IR pole. However, 
this is not to be trusted, since for positive bo we should interpret our results only as 
a power expansion in as, as discussed earlier. 

In case the behaviour of GikP') — G'(O) is of the type k'^lo^ik/Q) instead of a 
simple power, it is easy to convince oneself that the l/Q^ correction will be enhanced 
by n inverse powers of ag- In fact, the inclusion of powers of log(fc/Q) can be achieved 
from formula (|2.8| ) by taking derivatives with respect to p. Thus, since the large order 
behaviour of the expansion of Ip has the form 

n+l 



OC Tin + V 



( 2boac 
V P 



(2.11) 



-6- 



taking a derivative with respect to p we get the leading behaviour 

oc '-Tin + 2) ( ^V^' = -^Tin' + 1) ( ^V'^' (2.12) 

with n' = n + 1, which corresponds to a l/og enhancement. 

We have therefore found that in our approach the coefficient of the power correc- 
tion will depend upon the behaviour of S'(</), k"^) for small fc^, which is to say, upon 
how the definition of the shape variables for 4 partons goes to the 3-parton definition 
in the collinear limit. This indicates that the coefficient of the power correction can- 
not be simply factorized in terms of the three-body definition of the shape variable, 
a fact that we will examine in more details in the following sections. 

3. Details of the calculation 

We begin with some kinematical preliminaries. Outgoing legs for the three- and 




Figure 2: Labeling of external lines for three- and four-parton processes. 

four-parton process are given in fig. |^. We will call pi the momenta of the outgoing 
legs, Ei their energies, and the invariants will be defined as 



Q = V(^^0 

Sijk = {Pi+Pj+Pkf, yijk = ^- (3.1) 
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For the 4-parton process we have 



- ^234 ^ Q^^im. (3.2) 
In the three-body case these simphfy to 

The maximum value of S34 for fixed S134 and S234 is reached when pi and p2 are parallel 
and opposite: 

.34 = (3.4) 

Defining 

= = 1 - 2/234 , = = 1 - 2/134 (3.5) 

we have the constraint 

S34<g'(l-a;i)(l-a;2). (3.6) 

We follow here the calculation of ref. 0, from which many of the following results 
are taken. We work in ci = 4 — 2e dimensions. The three body cross section is given 
by the formula 

rfa(^) = H^C, ( ^^T(xi,X2) 
27r V Q2 y r(l -e) ^ ' ' 

9{xi + X2 — 1) ((1 — a;i)(l — X2){xi + X2 — 1))~^ dxi dx2. (3.7) 

The constant H is the normalization of the Born 2 body cross section, o"*-^-' = H. The 
four-parton cross section is 

, (A) TT^^sCp /47r/i^\'' 1 

^ ^ ^2Tr\Q^] f(rr^((^ ~^i''(^^^2) -2/34) (2/34 + + X2 - 1) 

^((l-xO(l-X2)-2/34) ^(2/34 + Xi+X2-l) rfx^ rfX2 "^^^ ^'"^^ ^ ^ ' 



27r V r(l-e) 

(^(1-^))-^ rf;^-i^sin-2^0W % T(^)(xi,X2, 2/34,^,^0 (3.8 

^^6» 2/34 

where 

T(^)(xi,X2,2/34,2;,6'') = T^t\\{^i,X2,z,e') + V{xi,X2,yM.z,e') (3.9) 
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1 - e ; 

R{X,,X2) = ,^ . ' r 3.12 

(1 - Xi)(l - X2) 

■^sin-^'e'de' = 1. (3.13) 

The V term (which can be extracted from ref. |^) vanishes for 1/34— ^0. We have 
kept the four-parton phase space factorized into a three-parton component, describ- 
ing the production of a quark, an antiquark and a gluon, and a two-parton term, 
corresponding to the decay of the virtual gluon into a quark-antiquark pair. 

In order to study the infrared renormalon, we must now include all the vacuum 
polarizations in the three and four-parton processes. First of all, we need a formula 
for the vacuum polarization in the MS scheme. We obtain 



U{k') = -i- 



37r \ /i^ / e 



r^ri-e) 1 



N(e) = U7rYT(l + e)^ — ^- f— ^ 1. (3.14) 

V ; V ^r(i_2e)(i_|e)(i_2e) ^ ' 

It is now easy to show that insertion of all vacuum blobs into a gluon line amounts 
to the following replacement 

as-^ . (3.15) 

The renormalizability of our theory implies that all divergences are removed by a 
redefinition 

as^Z ftren (3.16) 

where in the MS scheme Z is a power expansion in aj-e^, whose coefficients contain 
only inverse powers of e. After renormalization, our replacement rule will then become 

tts — ^ : ^ — 3.17 
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so that we must have 

= 1 + a,en&o -■ (3.18) 

Therefore, the resummation of all bubbles plus charge renormalization amounts to 
the replacement 

«s r—^. — (3-19) 



1 - asbo 



We can now immediately write down the result for the three-parton cross section 
including the effect of all bubble insertions. In this case, in fact, the momentum 
flowing through the gluon propagator is exactly zero. Indicating with a tilde the fully 
resummed cross section we get 

= L— p. (3.20) 

Observe that in spite of the renormalization procedure we carried out, poles in e do 
remain, and they should in fact be interpreted as infrared poles, that will ultimately 
cancel against analogous contributions in the four-parton cross section. 

The case of the four-parton cross section is more involved. In this case one should 
remember that in the q;| factor one power of as should be complex-conjugated. We 
then have 

1 1 



1-asbo {y e"- N{e) - 1) i l-asbo{y e"^^- N{e) - 1) i 



= da^"^^- ^ : (3.21) 



where we have defined 

y^(^y. (3.22) 

In order to make the infrared cancelation explicit, we will proceed as follows. Let us 
call G a generic infrared safe jet shape variable. For our purposes, G is characterized 
by two functions 

G(3)(xi, X2), G'(^) (xi, X2, y34, 0', z) (3.23) 
and infrared safety will imply that 

]im G^'^\x^,X2,ys4,0',z) = G^^\x^,X2). (3.24) 

2/34 — ^0 
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We are implicitly assuming that G does not receive contributions from the two-parton 
final state, and therefore it has a power expansion that starts at order a^, which is 
the case for all shape variables usually considered in e^e~ physics. The value of G in 
our model will be given by 

G = J da^^^G^^\xi,X2) + J dd'-^^ G'^^\xi,X2,y3A,0',z). (3.25) 

We now rewrite the above expression in the following form 

G = Gvirt + G4 + Gy + GcoU (3.26) 

with 



Gvirt = j dd^^^G^^\xi,X2) 
G4 = /rfa^") [G'^^\xi,X2,y3A,e\z)-F{xi,X2,yM)G'^''\xi,X2) 
Gv = Jda^^^ F{xi,X2,y34)G'^^\xi,X2) 

,3^2,2/34) G'^^^(xi,X2), (3.27) 



and 



N 1 - Xi 1 - X2) - yu) 1/34 + X1+X2- . 
F{Xl,X2,y34) = 77Z rrz 7777 7T^ e{Xl+X2-l) 

((1 - xi)(i - X2)r [xi + X2- ly 

rfa(^) = da'^^ + daiin- (3.28) 

The separation of da^^^ into the collinear and V term is performed according to 
eq. ( |3.9| ). The F factor is chosen in order to simplify the j/34 integral in the Gcou term. 
The G4 and Gy terms are free of infrared singularities, because the integrands vanish 



for 1/34— s>0, so that e can be safely replaced by in their expressions. In ref. [|I0[ it was 
shown that after 6' and z integration Gy is of order 7/34. In GcoU, the integration in 
z, 9' and 2/34 can be performed, because G*-^^ does not depend upon these quantities. 
The z and 9' integrals are easily done. The R term vanishes after angular integration, 
and the z integration gives 

,,(,(i^,))-.£!±(i^i)!^-r^(i-) 2(1-.) _ 2iV(.) 



r(l - 2e) (1 - 2e)(3 - 2e) 3 (47r)T(l + e) ' 

(3.29) 
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Using the identity 

we compute the yu integral 

r-{l-a;i)(l-X2) dy^^ 



dyzA dy ( Q"^ 



(3.30) 



(1 + 
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1 + 



asboN{e) (l + ^ 
■ arctan 



arctan ■ 



(l + ^ 



bo 



sm evr 



sm en 

(1-Xi){l-X2) 



1 












)' 



sm en (^i^^ 
agbon 



cosevr 



N{e)y 



1 + asbo (log 



iV'(O) 



en 



(3.31) 



where we have set explicitly e = in the first term. Using the identity 

en 



sine7rr(l + e)r(l - e) 



(3.32) 



we write our integral as 



Gcou = - [ -rr^d {xi + X2-1) dxi dx2 
J In 



X^T(xi,X2) G(=^)(xi,X2) 



37r 



arctan ■ 



OLsb^n 



(l-a'l)(l-X2)Q^ 



asCp I Anfi"^ 



2n \ J r(l-e 

c^ sTr 1 

3n as bo + 



1 + asbo log ■ 

6 (Xi + X2 - 1) [(1 - Xi)(l - X2){xi +X2- 1)]"' 

T{xi,X2) G^^\xi,X2) dxi dx2 (3.33) 



where we have defined according with ref. 11 



1 1 

— = _-iV'(0). 



(3.34) 



Remembering the identity for b^ we see that the second integral in the above equation 
cancels exactly Gvirt, so that all infrared divergences cancel, and we can write 



G — Ga + Gv ~\~ G<- 



Gs — Cvirt + Ccoll 



(3.35) 
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The Gs term is given by the first term of eq. (|3.33|) , and the 5* suffix stands for 



"Sudakov", since (in some sense) this term comes from the incomplete cancelation 
between the real and virtual diagrams, induced by the running of a^. It can be 
written as 

zn Jxi+X2<i Oq tt 

where 

kl = {l-xi){l-X2) Q\ (3.37) 



4. The Sudakov term 

A term of the form of Gs appeared first in ref. and was there used to param- 
eterize the 1/Q correction to shape variables like the average value of 1 — t, where t 
is the thrust. In ref. it was argued that the power corrections of the order 1/Q 
factorize in the form of eq. ( |3.36 ) for a generic shape variable, as well as for other 



processes. In our calculation, this term, before the Xi,X2 integration, does not have 
any renormalon, since it is an analytic function of ag near the origin. When we inte- 
grate over xi,X2, and approach the singular two-jet region, non-analytic behaviour 
may arise. This is due to the fact that we are integrating over the Xi, X2 values where 

l/as{kl) = l + asbo\og^-^0. (4.1) 

Our formula differs from the one of ref. P] only by the replacement 

(^sil"^)^-, — arctan(ag(g^)6o7r) . (4.2) 
Oovr 

In our case, however, there is really no singularity when we integrate over the Landau 
pole, since the arctangent is a bounded function. Let us compute the contribution of 
the Gs term to the average value of 1 — t. For a three massless body system, thrust 
is simply 

t = max(xi, a:2, 2 — xi — X2). (4.3) 

The leading contribution to 1 — t comes from the region where xi and X2 are very 
near 1, so 



t)s = ^2/ dx,dx2 ^^^-^'^ arctan(a,(fci)6ovr) 



2tt Jxi>x2 (1 - a;i)(l - 0:2) bon 

2Gf dy a.Tcta.n{ds{yQ^)bon) 



y/y boTl 



(4.4) 
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In order to evidentiate the structure of the infrared renormalon in the above formula, 
we compute the integral 

^ dy paYctan{as{yQ'^)bo7i) 
— ?/2 4.5 

y OqIT 

for arbitrary p > 0. We integrate by parts, and obtain 

^ 2 arctan(Q;5(Q)6o7r) ^2 dy e bo 



pboTT pJo y -|(^ + (M' 

2arctan(as(Q)6ovr) , 2bo . 
= + — Ip. (4.6) 

pboTT p 

The first term is analytic in as near the origin, while the second term, given by 
eq. ( |2.10 ) has an infrared renormalon located at z = p/{2bo), which corresponds to a 
power correction. For the case of (1 — t) we found therefore a 1/Q correction. 



5. The four— parton integral 



The terms G4 and Gy cannot easily be done analytically, because they depend in 
an intricate way on the four-parton phase space. Observe that 



G.+G 



V 



rfaW G'W(xi,a;2,2/34,^',^) -^(xi + a;2 - 1) rfa^i G^'Ha^i, ^2 



dy G4v{y) 



a. 



(l + as&olog^) +(asM' 



(5.1) 



where we have defined 

GMy) 



Ky- yu)[d(y^*'' G^'^\xi,X2,y3i,9',z) 
-e{xi+X2 - 1) da^^iiG^^\xi,X2) 



(5.2) 



It is clear that the small y behaviour of G4v{y) controls the power correction due to 
the IR renormalon. In particular, if 

,P/2 



GMy)'-^'^^- 



y 



(5.3) 



the position of the renormalon will be at p/ (26o), corresponding to a power correction 
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L .O t 0\ t } 




1 QQ7(^^ 

L.JJ t \o } 








9^ '^l (7^ 


1 DQ f>(A^ 






P 


* 


* 


51 Si's"! 


797l'4'l 


801 


849l'6'l 


mjj- > 0.1 


.384(2) 


-30.22(7) 


-62.8(3) 


-66.0(5) 


-13.60(15) 


-.035(10) 


P 


* 


* 


1.365(4) 


1.956(7) 


3.372(11) 


7.2(2) 


> 0.1 


.252(4) 


-119.6(3) 


-1856(5) 


-5250(30) 


-14610(140) 


-44600(800) 


P 


* 


* 


-.381(3) 


1.098(5) 


1.110(9) 


1.031(17) 


c > 0.2 


.934(3) 


-21.11(11) 


-110.5(10) 


-129(3) 


-145(11) 


-150(30) 


P 


* 


* 


.563(9) 


1.86(2) 


1.90(7) 


2.0(2) 


EECcut 


.4150(18) 


-2.12(2) 


-13.9(2) 


-51.8(13) 


-180(8) 


-630(40) 


P 


* 


* 


.371(16) 


.85(3) 


.92(4) 


.92(7) 



Table 1: Results for G^viv) for various shape variables, for y = 10^^, . . . 10^®. 
The line marked p, in the column corresponding to y = lO"'^, is the exponent one 
would obtain from the above table by fitting the pair of numbers on the line above, 
corresponding to y = lO"'^"'"^ and y = lO"'^ with the form y^^'^/y. 



We computed G^y numerically for ?/ = 10 ■'^ j = 1, . . . 6, for the following shape 
variables: (1— t), (6'(0.8— t)), {rrijj), where rajj is the heavy jet mass-squared according 
to the thrust definition, (^(m|^ - 0.1)), (^(o-O.l)), (^(c-0.2)), and for the weighted 
average of the energy-energy correlation away from the back-to-back region 

EEG^^t = EEG{cose) sin^^c/cos^. (5.4) 

J-0.5 



For the exact definition of these quantities, see ref. [|T^. The results are given in 
table m. For each value y = 10~^ , we also give the power p that is obtained by fitting 
G^viy) in the two points lO^^-'^^-' and 10"-' with a function proportional to y^^'^/y. 
As anticipated in the previous section, we can see from the table that the four-parton 
contribution can give 1/Q power suppressed corrections to quantities like the average 
value of 1 — t. 

It is easy to identify regions of integration that give such type of contributions. 
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Consider for example the region 

v/yil < (1 - xi) < VyiKi + ?7) Vyil < (1 - ^2) < Vy34(i + ^) (5.5) 

for 77 <^ 1 and independent of 2/34. In this configuration, 1 — i is always of order ^/yal. 
For small ^34 the emitted gluon is soft, so that the amplitude factorizes in terms of the 
three body amplitude in the soft limit, and a term depending upon the orientation of 
partons 3 and 4. The three body amplitude gives an integral of the form 

(5.6) 

JVy^ 1- xi 1- X2 

independent of ^34. Next, we have to weight the amplitude with 1 — t, which is of 
order ^/yu, and integrate in dy^i/yzi- This corresponds to G^v oc y/yu/y34., which, 
as we have seen, yields a 1/Q correction. 

One may wonder whether one may still recover a factorized form, similar to the 
term Gs, by a suitable redefinition of the effective coupling. Looking at the quantity 
{mjj) we see that this is not the case. In fact, the invariant mass of the heavy jet 
is equivalent at the 3-parton level to thrust, namely we have 1 —t — mjj/Q'^. This 
identity is no longer valid at the 4-parton level. For example, in the soft configuration 
we have just considered, where partons 1 and 2 are back-to-back, and so are partons 
3 and 4, it is a simple exercise to show that the relation is instead 1 — t — 2mjj/Q'^. 
Therefore, any three-body factorization formula would fail in this simple case. Shape 
variables that are identical at the 3-parton level, but differ at the 4-parton level, have 
different coefficients for the leading power correction. 

Let us now focus on shape variables that depend upon final state configurations 
that are far from the two jet region, for which the Sudakov term does not provide a 
leading 1/Q power correction. For some of these variables, e.g. (^(0.8 — t)), {0{mjj — 
0.1)), {9{c — 0.2)), we see no evidence for power corrections of the form 1/Q, but 
instead we find a l/Q^ correction (in the case of {0{mjj — 0.1)), the term has a 
rather small coefficient, so this behaviour does not become apparent until y ^ 10~^, 
where G^vin) becomes constant at +1.43, indicating p — 2). For {9{o — 0.1)) we 
instead observe a. 1/Q type of correction. A leading 1/Q correction is also observed 
for EECcnt, in spite of the cut that avoids the back-to-back region. As wc will see 
shortly, this is due to the fact that EEC receives contributions from configurations 
near the two jet region also for angles far away from the back-to-back configuration. 

These findings can be easily justified by examining the singular integration region 
for various quantities. First of all we will consider thrust. Let us look at parton 
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configurations near the coUinear limit. By kinematical reasoning one can convince 
oneself that the thrust axis is either along parton 1, parton 2, or along the sum of 
partons 3 and 4, and the thrust is given by xi, X2, or (2 — Xi — X2) — '^^y-a/t + 0{y 



respectively. In all cases, it differs from the thrust of the corresponding configuration 
with 1/34 = by terms of the order of 1/34 or less. This behaviour gives rise to a 1/Q^ 
power correction. In the case of oblateness, we can instead identify a region where 
a ^2/34 behaviour arises, leading to a 1/Q power correction. One such configuration 
is depicted in Figure H. One projects the event onto the plane orthogonal to the 




a. 



o, 




thrust projection 



Figure 3: Oblateness in the collinear limit. 

thrust axis, and then oblateness is defined as the difference between the major and 
minor axis. For the particular configuration shown in the figure, oblateness is just 
the difference between the distance 1-2 (called /-major) and the distance 3-4 (called 
/-minor) in the projected event. It is easy to convince oneself that the 3-4 distance 
is proportional to ^Jy^. In the limit of 1/34— >0 oblateness behaves therefore hke y/j/il, 
which generates a, 1/Q correction. 

The 1/Q correction in the case of the energy-energy correlation has instead a 
very different origin. One such correction arises from the Sudakov term. In the 3- 
body configuration, when parton 3 is soft, the Sudakov term for the energy-energy 
correlation is 

Cp 1 pQ^i'i^ dk±_ arctan(as(/ci)6o7r) 



sai^e EEC {cos e) = 8- 



(5.7) 



2n sin 9 Jo Q bon 

The integral over k± generates a 1/Q correction. Other contributions come from the 
G4V term. Consider some weighted average of the energy-energy correlation, with 
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a weight f{6) in an angular interval that does not include the back-to-back region. 
Then, for example, the contributions coming from partons 1 and 3 would be 

J da,f{e)E,E,. (5.8) 

If parton 3 (the gluon) splits into a quark-antiquark pair, carrying fractions z and 
1 — z oi parton 3's momentum, and having an opening angle uj. Also assume for 
simplicity that the splitting takes place in the 1-3 plane. The contribution in the 
coUinear limit is then 

j da, El ^3 {f{e + zuj){l -z) + fie - (1 - z)uj)z) 
^ J da, E, E, [f{6) + i f"{e)uj' z{l - z)^ 

= J da, E, E, ^f{e) + ^ ne)su/El^ , (5.9) 

which seems to give rise to a l/Q^ power correction. This is in fact not the case, since 
we integrate over the region where E^ becomes small. In this region the cross section 
behaves as dE^/E^, so the E3, integral in eq. ( |5.9| ) yields 




(5.10) 



and we see that in this way al/Q power correction does arise. Therefore, the energy- 
energy correlation receives 1/Q power corrections for all values of the angle, coming 
from the region near the two-jet limit. It might appear that if we set f{6) = constant, 
this 1/Q correction would be zero. However there would still remain ^-function 
weights defining the edges of the integration region, which would give equivalent 
1/Q terms. 

Unlike the case of oblateness, we see that for the EEC the correction arises because 
the kinematic region near the two-jet configuration contributes for all values of the 
angle. All other commonly considered shape variables depend instead upon the three 
jet region for intermediate values of the shape parameter. 



6. Higher order terms 



We have seen from the previous section that, except for special cases like oblate- 
ness, 1/Q corrections arise from configurations with a soft gluon emission, where the 
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gluon virtuality is of the order of its energy, followed by the decay of the virtual 
gluon into massless partons. This process cannot occur away from the two jet region 
at leading nf, but it can certainly arise at subleading Uf. We may imagine adding 
to a three-jet, qqg configuration a soft, off-shell gluon (i.e. with energy of the same 
order as its virtuality) decaying into a massless parton pair. It is difficult to imagine 
any shape variable that will not receive 1/Q corrections from this kind of process. In 
fact, shape variables are typically linear in the parton momenta, as dictated by the 
requirement of insensitivity to coUinear splitting. The production of a soft, off-shell 
gluon reduces linearly the energy available to all the other partons, which in general 
may affect the shape variable linearly in the gluon energy. Since the cross section 
for soft gluon emission has the characteristic behaviour dEg/Eg, and the emission 
coupling will be evaluated at the virtuality of the gluon (assumed to be of the same 
order as Eg) it follows that 1/Q corrections are present. We have not, of course, 
rigorously proven this fact. Needless to say, if shape variables that never develop 
1/Q corrections were found, their importance for the determination of as would be 
enormous. 

Let us therefore assume, for a moment, the pessimistic (and perhaps realistic) view 
that shape variables always develop a. 1/Q correction at some order in perturbation 
theory. Let us consider, for example, thrust with a cut t < 0.8, so that we are always 
in the three jet region. According to the above argument an extra soft gluon emission 
will generate a correction of order 1/Q. It seems plausible however that the hard 
real emission contributes a factor of as{Q'^), such that the overall correction is of the 
order of as{Q)/Q. This would again be a very important fact. It would tell us that 
some shape variables are indeed better than others, in the sense that their 1/Q power 
correction carries an extra as{Q) suppression. 

It may also be possible that the 1/Q suppression will turn out to be enhanced 
by a power of \og{Q/A), which would compensate the as suppression. This could be 
produced, for example, by a 5-parton term that behaved as y^yis logy45 when particles 
4 and 5 become coUinear. Whether these logarithmic enhancements are present or 
not is a matter that ought to be clarified with further studies. In the present work, 
we simply remark that it is conceivable that one may find shape variables in which 
the enhancement is not present, and that therefore do have a as suppression of the 
1/Q power corrections. 
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7. Conclusions 

In the present work, we have proven that even in the simple model of QCD at large 
rij, shape variables in e^e~ annihilation show remarkably different properties with 
regard to power corrections originating from infrared renormalons. In particular, we 
have shown that in the large nj limit, variables like {1 — t), {rri^), and the EEC for 
any value of the angle, develop a. 1/Q correction, while thrust and the c parameter do 
not develop any 1/Q correction in the region where the two jet configuration does not 
contribute. Another remarkable result is that oblateness develops a 1/Q correction 
even away from the two jet region. 

We compare our findings with the results of refs. and . We recover a correc- 
tion term with the factorized form proposed there, but we also find an extra correction 
that spoils factorization, since it specifically depends upon the 4-parton definition of 
the shape variable. Thus, two shape variables that have the same 3-body expression, 
but differ at the 4-body level, will have different power corrections. The discrepancy 
with the authors of ref. can be tracked back to the fact that they assume to some 
extent the validity of the perturbative expansion even when the scale of ag is very 
low, a fact that does not take place in our calculation in the leading Uf limit. 

We argue that even shape variables that do not develop a 1/Q correction at the 
leading large Uf level, may develop one at subleading level, and therefore in the full 
QCD. We conjecture that the leading power correction to shape variables will have 
in general the form ag((5)/(5, and one may classify shape variables according to the 
value of n. It may therefore be possible to find a class of shape variables with leading 
power correction of the form ag{Q)/Q. With these shape variables, the influence of 
non-perturbative effects upon the determination of would be truly negligible at 
LEP energies. 
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